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ABSTRACT 

We present SOAR/OSIRIS cross-dispersed near-infrared (NIR) integrated spectra of 12 
Galactic globular clusters that are employed to test Maraston (2005, hereafter M05) NIR 
Evolutionary Population Synthesis (EPS) models, and to provide spectral observational con- 
straints to calibrate future models. We measured Equivalent Widths (W\) of the most promi- 
nent NIR absorption features: A 1.49/xm, Mgi A 1.58/xm,Fel/Mgi, A 1.59^m, Sii, A 1.71/im, 
Mgl, A 2.21/im, Nal and A 2.26/mi, Cal as well as the A 1.62/im, A 2.29//m, CO and A 
2.05 A*m, CN molecular bands. Optical W\ of G-band (4300A), H/3, Mg 2 , Fe I (4531 A, 527A 
and 5335A), and Na I (5897 A) were also measured. The globular clusters W\ were compared 
with model predictions with ages within 4 — 15 Gyr, and metallicities between ^ ZQ and 
2 ZQ. Observed integrated colours (B — V , V — I and V — K s ) were also compared with 
models. The NIR integrated spectra among our sample appear qualitatively similar in most the 
absorption features. The M05 models can properly predict the optical W\ observed in globu- 
lar clusters. Regarding the NIR, they do underestimate the strength of Mg I 1.49/im, but they 
can reproduce the observed W\ of Fe I 1.58/xm, Si I 1.59/im, and CO 2.29 /im, in about half 
of our sample. The remaining objects require the inclusion of intermediate-age populations. 
Thus, we suggest that the presence of C- and O-rich stars in models is important to reproduce 
the observed strengths of metallic lines. Another possibility is the lack of a-enhancement in 
the models. In the case of the optical and NIR Fe I lines, standard models and those that in- 
clude blue horizontal branch stars, produce similar results. A similar trend is observed for Na I 
5895A, while in the case of the G-band, the models with blue horizontal branch do describe 
better the observations. For most of the sample the optical to NIR colours are well described 
by the M05 models. In general, M05 models can provide reliable information on the NIR stel- 
lar population of galaxies, but only when W\ and colours are taken together, in other words, 
W\ and continuum fluxes should be simultaneously fitted. However, the results should be 
taken with caution, since the models tend to predict results biased towards young ages. 

Key words: Globular clusters: general, Galaxy: stellar content, Near Infrared 



1 INTRODUCTION 

The study of stellar populations is a critical step to understand 
the continuum emission of galaxies, even in active galactic nu- 
clei (AGN). This occurs because several components such as the 
non-thermal continuum, dust emission, the presence of an active 
nucleus, and the stellar population itself, sum up to the integrated 
spectrum. By analysing the stellar content, information can be ob- 
tained on critical processes like recent star formation episodes, the 
evolutionary history of the galaxy, and the connection between the 
active galactic nucleus and starburst activity. 
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Observations in the near-infrared (NIR) are important to study 
the integrated stellar populations of galaxies, since this is the most 
convenient spectral region accessible to ground-based telescopes 
to probe highly obscured sources. However, tracking the s tar for- 
mation in the NIR is complicated ( Origlia & Oliva 2000). Nev- 
ertheless, NIR stellar absorption features are widely believed t o 
provide means for recognising red supergiants Joiiva et alii 19951) . 
since they arise as prime indicators for tracing starbursts in galax- 
ies. Besides the short-lived red supergiants, the NIR also includes 
the contribution of thermally- pulsating asymptotic giant branch 
(TP-AGB) stars, enhanced in young to intermediate age stellar pop- 
ulations (0. 2 < t sC 2 Gvr, lMarastoiJl998l.l2005l : lRiffel et al.l200l 
l2008ll2009h . The TP-AGB phase becomes fully developed in stars 
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with a degenerate carbon oxygen core (see llben & Renzinlll 19831 
for a review). Evidence of this population in the optical is usu- 
ally missed, as the most prominent spectr al features assoc iated with 
this evolutionary phase falls in the NIR <Marastorjl2005l . hereafter 
M05). 

Usually, the interpretation of a galaxy's stellar population 
involves a synthesis approach that, in general, is based on 
the mixing of simple stellar populations of different ages and 
metallicities that provides the best match to the galaxy spec- 
tra. The combinations can be done usi ng equivalent widths (W\) 
and selected continuum fluxes (e.g. iBica & Alloinl [19861 : iBical 
1 19881 ; Riffel et al. 200 8]) or by fitting the wh ol e underlying spec- 
trum (e.g. Ibid Fernandes et al] 12004 12005 jlbt lAsari et al ]|2007l: 
ICid Fernandes et alj 120081 : iRiffel et ail 120091) However, in both 
methods the most important ingredient in the SP synthesis is the 
spectral base. An ideal base should cover the full range of spectral 
properties that occurs in a galaxy sample, providing enough resolu- 
tion in age an d metallicity to properly address the desired scien tific 
question (see ISchmidt et al1ll99ll : ICid Fernandes et al .11200511 for 
example). 

A reliable base to probe the stellar population of galaxies is 
a library of integrated spectra of star clusters such as the one con- 
structed by IBica & Alloinl £[986). However this base is restricted 
to the optical region and metallicities lower than solar ([Z/Zq] < 
0.1). The main advantage of this approach over those based on 
spectra of individual star s, the so called Evolutionary Population 
Synth esis (EPS) models dBruzual & Charlotll2003l : IVazdekis et al] 
I2010L M05), is the reduced number of variables. While the latter 
are essentially described by the temperature, gravity and metallic- 
ity, u ncertainties in the form er are reduced to age and metallicity. 
Also, IBica & Alloinl {T986) method is free from any assumptions 
on stellar evolution and the initial mass function. However, one ad- 
vantage of the EPS models is their large range in ages and metal- 
licities, while observed spect ral libraries are restricted to the opti- 
cal (e.g. lBica& Alloinlll986l) or cover only small fractions of the 
NIR spectr al region, like the region a round the 1.6^m and 2.29/xm 
CO bands (lOriglia et all 1 1 994lTl997h or only the if -band spectra 
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The above considerations justify a project to check the reli- 
ability of NIR EPS models available in the literature and to set 
important constraints on the absorption features observed in this 
spectral region. We focus here on the observation of NIR (~1.2/iin- 
2.35/im) integrated spectra of globular clusters. By their very na- 
ture, globular clusters spectra should be reproduced by simple stel- 
lar population theoretical models. 

This paper is structured as follows: In section [2] observations 
and data reduction procedures are described. Results are presented 
and discussed in section[3] and conclusions are presented in section 
II 



2 SAMPLE SELECTION, OBSERVATIONS AND DATA 
REDUCTION 

As part of an ongoing project to investigate NIR spectral properties 
of globular clusters, w e have selected a repr esentative sub-sample 
of the star clusters from lBica & Alloinl (fl98a) to test NIR EPS mod- 
els. Our original sample was composed by 27 star clusters covering 
almost all the possibilities of the age/metallicity that occur in galax- 
ies. 

As stated by M05, even bright optical star clusters have the 
NIR light dominated by stars in short-lived evolutionary phases like 
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the bright red giant branch (RGB), as well as the TP-AGB. Thus, to 
avoid the resolution of individual stars, two strategies were used: (i) 
the telescope was defocused and (ii) non sidereal guiding, sweep- 
ing the area corresponding to the cluster centre. Consequently, inte- 
grated spectra of the mixed light of a representative fraction of the 
cluster stars were obtained. 

Using these strategies, cross-dispersed (XD) NIR spectra of 
23 star clusters were obtained from October 2006 to June 2009, 
at the Ohio State Infrared Imager/Spectrometer (OSIRIS) attached 
to the 4.1 m Southern Astrophysical Research (SOAR) Telescope. 
The detector is a 1024 x 1024 HgCdTe array. At the XD mode 
the spectrograph provides simultaneous coverage of the region be- 
tween 1.2jitm and 2.35/im with a resolving power of R~1200. It 
is interesting to note that in this mode it is possible to observe the 
1.2/im to 2.35/^m region in a single shot, avoiding the seeing and 
aperture effects when observing in single bands. 

Because of the extended nature of the science targets, it was 
necessary to take separate sky exposures, following the pattern 
ABAB, where A is the on-source observation, with the object cen- 
tred on the slit, and B represents off-source exposures. To remove 
telluric lines and to flux calibrate the spectra, standard stars were 
observed, following the nod pattern along the slit. 

Given the nature of the NIR light emitted by globular clus- 
ters, the above observational strategy and configuration produced 
only 12 reliable (S/N > 25) integrated spectra of globular clusters 
(TabQJ. We list their observed colours and some basic properties 
in Tab . [2] As it can be seen from Tab. [2] and Fig. 1 of IBica et al] 
(2006), seven clusters are metal-poor, while the metal-rich ones are 
almost all located in the Galactic bulge (dGC ^ 5kpc). 

The spectral reduction, extraction and wavelength calibra- 
tion procedures were performed using XD-SpresQ IRAF software 
task (Ruschel-Dutra et al. 2010, in preparation). This tool fol- 
lows the standard NIR cross dispersed reduction procedures, as 
listed: (i) to remove sky emission lines from the integrated spectra, 
sky exposures were subtracted from the object, (ii) the resulting 
images were summed up and then divided by a normalised flat- 
field image, (iii) extraction was performed following the standard 
echelle procedures, (iv) wavelength calibration was based on the 
OH lines present in th e sky exposures with the values given by 
lOliva & Origlial jl992). (v) the 1-D spectra were then corrected for 
telluric absorption by comparison with an A0V star spectrum. The 
stellar atmospheric absorption lines were identified and removed. 
This procedure was done using the TELLURIC task, (vi) finally, the 
flux calibration was achieved by fitting a black-body curve to the 
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Table 2. Observed colours and some basic properties of the cluster sample. 
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Table Notes: (a) From Bica et aljj2006h : (b) From - http://www.physics.mcmaster.ca/Globular.html I Harris 1996); (c) From lCohen et alj 12007). 
The NIR colours are reddening-corrected ICohen et al.l2007h while the optical ones are uncorrected faarrij 19961) . 



standard star featureless spectrum, using the task CALIBRATE of 
the IRAF software. Final reduced spectra, are shown in Fig.Q] 



3 RESULTS AND DISCUSSION 
3.1 General overview of the spectra 

Spectra of the 12 Galactic globular clusters covering simultane- 
ously the NIR spectral region between 1.2/xm - 2.35/xm are pre- 
sented here for the first time. The full NIR spectral energy distribu- 
tion (SED) of the clusters is shown in Fig.Q] The NIR SED of the 
clusters is rather similar among the 12 objects analysed, decreasing 
smoothly in flux with wavelength, and present similar absorption 
features. The shaded area in Fig.Q]denotes a region highly affected 
by telluric absorption, thus it was left out of further analysis. 

In order to highlight the absorption lines, the spectra were lin- 
earised by their continuum. The result of this process is shown in 
Figs.|2]to|4] Clearly, the spectra of the globular cluster sample ap- 
pear qualitatively similar in most of the NIR absorption features. 
Such similarity can be observed in Figs.[2]to[4] where many atomic 
absorption features like A 1.49/im Mgl, A 1.58/xm Fe I/Mgl, A 
1.59^m Si I, A 1.71^mMg I, A 2.21/xm Nal and A 2.26/xm Cal as 
well as the A 1.62pim, A 2.29/xm CO and A 2.05/im CN molecular 
bands are clearly detected and identified in the spectra. The centre 
and width of the strongest lines are identified in Figs. [2] to [4] The 
similarity between the spectra is not surprising, since the sample is 
composed only by globular clusters and thus, the integrated spectra 
is dominated by stars of relatively similar spectral type. 

It is also clear in Figs.|2]to|4]that the spectra show many weak 
absorption lines. However, as the S/N (see Tab.Q} is relatively poor, 
no effort was made to identify such weak absorptions. In addi- 
tion, the H and Jf-band lines pres ently detected are common in 
the brightest stars of star clusters l Origlia, Rich & Castro 2002; 
lOriglia et alj 12006) : iFrogel et alj|200ll ; IStephens & FrogelH2004l) . 
which indicates that our data reduction is consistent. 

To date the only effort in o btaining NIR integr a ted sp ectra of 
Galactic star clusters was from IOrigliaetal.l ( ll994l.ll997l). These 
autho rs used the IRSPEC infrared spectrometer (Moorwood et al. 
1 19911) attached to the European Southern Observatory (ESO) New 
Technology Telescope (NTT). Their instrumental setup allowed 



them to obtain long-slit spectra centred at 1.59^im (Si I), 1.62/iin 
(CO 6-3), and 2.29/xm (CO 2-0). These absorptions are clearly 
detected in our spectra (Figs. [2] - 14). For a compa r ison between 
both data sets, the W\ reported by lOriglia et ail i ll 994 1 19971) 
for the objects in common with our sample were plotted against 
those measured by us (Fig. [5}- The Pearson correlation coeffi- 
cient (r 2 = 0.87) and a nonlinear regression were computed us- 
ing the FITEXY routine from the IDL Astronomy User's Library, 
which takes errors in both coordinates into account. We find Wa = 
-(3.66±1.65) + (1.93±0.29) W A ,i it , where W x means this work 
and W\,ut the literature data. A generalised correlation coefficient 
(CC=0.75), defined as the square root of the fraction of total y vari- 
ance explainable by the fitted function was also compu ted. Clearly, 
there i s a discrepancy between our values and those of Orig lia et al.1 
dl997l) . Probably, the difference in W x of Si I 1.59/wn and CO 
1.62/im occurs because these authors use rectified spectra, instead 
of locally defined continuum regions. For the 2.29/im CO band, the 
differences in W\ are associated with different line definitions (see 
Sec. |3.2t . I n addi t ion, the dif ference may be partly due to the fact 
that lOriglia et alj d 19941 1 19971) only sampled a very small fraction 
of the cluster core (6" x 4"). Thus, the random contribution of a 
few bright, red giant stars in such a small field, may strongly affect 
features in the integrated spectra. 



3.2 Equivalent Widths 

A straightforward approach for comparing a set of empirical spec- 
tra with models is by means of the W\ of their absorption features. 
Such a procedure is virtually free from flux calibration issues and 
reddening corrections, and may provide constraints on model pre- 
dictions. 

When dealing with optical absorption lines, the spectral win- 
dow of a feature is usually defined together with two continuum 
regions at the red and blue sides, which are used to trace a local 
continuum t hrough a linear fit to the mean values o f both continuum 
regions (e.g. lFaber et a l. 1985; Bi ca & Alloinl 19861) . Si milar defini- 
tions were made for som e NIR absorption features (e.g. lOliva et al.1 
ll995l : lFrogel et alfcOOll) , however, they were based in stellar spec- 
tra, where absorption features are typically a few kms -1 wide. As 
we are dealing with velocity-dispersion sustained (a v >10km s _1 ) 
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Figure 1. Globular Clusters NIR spectral energy distribution. The spectra were normalised by the mean value of their fluxes between 1.76^m and 1.77/imand 
shifted by a constant for display purposes. The shaded region was left out of the analysis (see text). The gaps represent regions strongly dominated by telluric 
absorptions, and/or coincide with band limits. 



clusters, and we aim to use these absorption features to investi- 
gate unresolved stellar populations in galaxies, we are forced to use 
spectral indices with broader band passes than those used in indi- 
vidual stars. Our continuum and band-pass definitions are listed in 
Tab. [3] They were taken from the literature, except for Mg I 1 .49 /im 
and CN 2.05^im, which we define in Tab. [3] In some cases it was 
also necessary to re-define the continuum regions in to avoid over- 
lap (Tab. |3j- However, special care was taken to use only regions 
free from emission/absorption lines. An example of these contin- 
uum and line limits is shown in Fig. [6] 

With the definitions for continuum and band pass listed in 
Tab. [3] W\ was measured for nine absorption lines, using the code 
PACCE, which computes W\ with user d efinitions for the line and 
continuum band-passes dVale et alj|2007h . The measured W\ val- 
ues are presented in Tab. [4] Errors w ere estimated according to 
Eq. 7 of lVollmann & Eversbergl feOOfj) . 

In addition, we measured on the 1 2 globular clusters, W \ us- 
ing the Lick indices definitions (e.g. IWorthev et al.l [1994J, taken 
from: http://astro.wsu.edu/worthey/html/system.htmll for the fol- 
lowing optical features: G-band (4300A), Fe I (4531 A, 5270A and 
5335A), H,S, MgH (5176A), and Nal (5895A), MgH (5176A). The 
W\ corresponding to the MgH absorption is called Mg2, which is 



Table 3. Line limits and continuum band-passes. 
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Table Notes: Table Notes: CBPb an d CBP r are the blu e and red band- 
pass boundaries for the conti nuum, (a) fromlRiffel et alj j2008h . (b) from 
lOriglia et all 1 19931) . (c) from lCesetti et alj J2009l) . 



adopted in the present p aper. The optical int egrated spectra of the 
clusters were taken from I Santos et alj d2002l) . Observation and re- 
duction procedur es, as well as data q uality of the optical spectra, 
are discussed in iBica & Alloir] Jl986h . The optical W x are listed in 
Tab. [4] Errors were estimated similarly to the NIR W\. 
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Figure 2. J-band final reduced spectra. The most prominent lines are Figure 3. H-band final reduced spectra. The most prominent lines are 
identified. The shaded region represents the line limits and the vertical identified. The shaded region represents the line limits and the vertical 
line shows the line centre. line shows the line centre. 



25 



3.5 



I 2.5 



1.5 



1.95 




2.1 2.15 2.2 
Wavelength (urn) 



Figure 4. K-band final reduced spectra. The most prominent lines are 
identified. The shaded region represents the line limits and the vertical 
line shows the line centre. 



W x =-(3.66± 1.65) + (1.93±0.29)W 14J , 
r 2 0.87 te CC 0.76 




"wi. (A) 

Figure 5. Comparison of the measured W\ with literature data from 
lOriglia et ail ll997t) . Features used are 1.59^m (Si I), 1.62/xm (CO 6- 
3), 2.29/im (CO 2-0). The Pearson coefficient (r 2 ), a generalised cor- 
relation coefficient (CC, see text) as well as the regression equation are 
shown. 
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Table 4. Measured W\ (in A) for the cluster sample. 
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104 3.50±0.24 1.53±0.36 1.33±0.18 5.26±0.04 1.09±0.16 U3±0.50 I.73±0.I6 1.59±0.31 

362 2.68±0.23 1.34±0.36 I.56±0.I8 1.33±0.04 0.72±0.16 0.83±0.51 0.63±0.I6 1.05±0.32 

1851 1.50±0.23 0.85±0.37 I.57±0.I8 1.38±0.04 1.09±0.16 0.52±0.5I 0.62±0.I6 173±0.31 

2808 2.46±0.23 2.00±0.36 I.83±0.I8 0.96±0.04 0.98±0.16 0.60±0.5I I.47±0.I6 1.15±0.34 

6388 2.70±0.23 2.14±0.36 1.3 1 ±0. 18 4.91±0.04 1.23±0.17 1.18±0.50 3.43±0.15 1.79±0.31 

6440 3.1I±0.25 2.24±0.36 1.8 1 ±0. 17 6.86±0.07 0.90±0.17 0.91±0.51 3.68±0.15 0.13±0.32 

6517 2.62 ±0.23 3.39±0.38 1.79±0.I8 — 0.97±0.16 0.85±0.51 2.42±0.15 0.92±0.32 

6528 3.04±0.28 3.04±0.36 0.58±0.I8 8.07±0.05 1.7I±0.17 1.72±0.50 4.39±0.15 1.36±0.30 

6541 i.96±0.23 — 2.09±0.I7 0.94±0.04 0.99±0.16 0.89±0.50 I.46±0.16 0.80±0.31 

6553 9.95±0.33 1.94±0.45 1.71±0.19 8.27±0.06 2.71±0.17 1.73±0.50 3.50±0.16 — 

6864 1.26±0.24 0.73±0.37 2.49±0.I7 2.56±0.04 0.87±0.16 0.50±0.5I 0.78±0.I6 0.72±0.30 

7078 — 0.21±0.37 1.94±0.18 0.57±0.04 0.10±0.17 0.10±0.51 1.49±0.I6 — 



1.41±0.25 0.56±0.09 1.34±0.18 0.58±1.05 — 7.05 ±3.48 — 18.94±2.04 

5.98±0.24 1.38±0.09 0.74±0.22 1.97±1.22 7.10±0.16 — — 23.32±1.89 

0.64±0.24 1.08±0.08 — — — 0.52±4.61 — 4.46±2.07 

2.30±0.24 0.78±0.09 0.61±0.24 — 0.25±0.15 4.18±4.19 0.10±0.93 14.42±1.81 

4.69±0.23 >0.11 0.96±0.21 0.62±1.05 2.93±0.15 — 0.60±0.96 14.3 1 ± 1 .93 

2.34±0.22 0.42±0.08 I.26±0.I7 1.07±1.05 0.58±0.12 1.59±3.83 0.52±1.00 1 1.41 ± 1 .96 

3.22±0.22 1.36±0.08 0.69±0.18 — 0.59±0.13 — — 2.12±1.87 

5.14±0.23 2.67 ±0.08 2.89±0.21 1.93±1.03 — — 2.32±0.87 14.78±1.76 

4.07±0.24 1.74±0.09 — — — 0.64±3.62 — 9.85±1.86 

5.93±0.22 2.88±0.08 1.72±0.18 1.70±1.04 4.85±0.13 — 2.09±0.90 16.99±1.79 

5.00±0.22 1.71±0.08 — 0.78±1.06 0.64±0.19 — 2.64±0.90 3.76±1.95 

1.60±0.22 — — — 6.24±0.13 — — 2.3I±1.93 
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Figure 6. Examples of line limits and continuum band-passes. The labels 
(a), (b), (c), (d), (e) and (f) indicate the cluster from which they were taken, 
being: NGC 104, NGC 6388, NGC 6528, NGC 6440, NGC 7078 and NGC 
6553, respectively. 



3.3 Observations compared to models 

Theoretical spectral models for SSPs have become fundamen- 
tal in studies of the stellar population of galaxies. Thus, 
as more such models are available, it becomes increasingly 
important to test them over all the spectral wavelengths 
dCid Fernandes & Gon zalez-Delgado 2010). Our main goal in this 
paper is to compare the integrated spectra of actual globular clusters 
with theoretical predictions of M05 models. These models were 
chosen due to the fact that up to date they are the only ones that 
make predictions about the presence of almost all observed absorp- 
tion lines and bands in the NIR spectral region. 

In order to compare observations with models we measured 
in M05 models the same set of Wx (using our continuum and line 
definitions) as in the star clusters. The Wx of the optical absorption 
features were also measured on the models using the definitions 
described in Sec. 13.21 



Indices, such as Wx of spectral absorption lines, can be seen as 
a compres sed, but highly inform ative, representation of the whole 
spectrum dCid Fernandes 1 12007). Thus, they are one of the more 
suitable ways to compare the observables with model predictions. 
To compare observations with models, NIR spectral indices (this 
work), were plotted against each other and versus optical Wx- 
Some of the most significant correlations are shown in Fig. [7] to- 
gether with model predictions. 

In Fig.|7^ we compare the relation between NIR Fel 1.58 /im 
and the optical Fe I 4531 A with model predictions. The shaded re- 
gion represents the model values for ages from 4 up to 15 Gyr (dot- 
ted and solid lines represent the limits), with metallicities between 
Zq and 2 Z@ . Clearly, the models can reproduce the Wx of 
the optical Fe I absorption line of the star clusters (i.e. old ages and 
Z Zq). However, they only predict the values of the NIR Fe I 
line in about half of our sample. Similar conclusions are obtained 
when considering Figs.|7j5 and[7J;, in the latter case, however, for 
the NIR Si I line. 

ICid Fernandes & Gonzalez-Del gadol d201 Ol) have shown the 
importance of modelling the horizontal branch to properly estimate 
the origin of blue stellar populations in integrated light analyses. As 
M05 models consider the blue horizontal branch, we include such 
models in Fig. |7Ji (open triangles joined by a dashed line). Note 
that we have only included the 15 Gyr old population because the 
4 Gyr models are not available for all metallicities. No significant 
differences are observed in the prediction of blue horizontal branch 
and standard models in the case of the Fe I and Si I lines (Figs. [7^ 
toS). 

Similarly to Fig. [TJi, we compare the optical G band and the 
2.29 fim CO band in Fig. |7Jl. As in the case of the Fe I lines, the 
observed optical band values are almost all in agreement with the 
model predictions. However, in the case of the 2.29^im CO band, 
the models predict lower values than the observations for ~ 50% of 
the sample. A similar trend occurs for H/3 (Fig.|7j;) and Na I 5895A 
(Fig-Ef) . 

Regarding the presence of blue horizontal branch stars in the 
models, they do describe better the optical absorption, i.e. they pre- 
dict lower values for the G band. As expected, they also play an 
important role in determining the Wx of H/3. In contrast, the blue 
horizontal branch models follow the standard ones concerning the 
Nalline (Fig.|7f). 

The NIR Mg I line and the optical Mg2 are compared with the 
models in Fig.|7g. As for the Iron lines, the optical feature is prop- 
erly predicted by the models, while they fail in the case of the NIR 
Mg I line, since in this case the predicted values are about one or- 
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Figure 7. Comparison of the measured W\ with model predictions. The shaded region represents models predictions for ages from 4 up to 15 Gyr (dotted and 
solid lines represent the limits), with metallicities between — i-jj ZQ and 2 ZQ. Open triangles joined by a dashed line represent models which consider the 
blue horizontal branch. Asterisks joined by the dot-dashed line represent models with 1 Gyr old. Metallicities are labelled, except the lowest value, which we 
left out for display purposes. 
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Figure 8. Same as Fig.|7]but for optical W\ and without the 1 Gyr models. 




Figure 9. Dependence of selected optical indices on [Fe/H]. Full and dashed lines are 4 Gyr and 13 Gyr old models, respectively. 



der of magnitude lower than the observations. The blue horizontal 
branch models follow the standard ones. 

In Fig. [7ji and |7J we compare NIR absorption lines against 
each other. It is clear that the models do not reproduce the observed 
strengths for Mgl 1.49/im. However, as discussed above they do 
describe the measured values for Fel 1.58/im and Si I 1.59/im in 
about half of the sample. As expected, the blue horizontal branch 
models are very similar to the standard ones. 

We also check the M05 models in the optical by comparing the 
W\ of optical absorption lines against each other in Fig.[8] Clearly, 
the models reproduce the absorption line strengths measured in al- 
most all objects. Although a dependence on metallicity is already 
apparent in Figs. 1817 1 we show this more explicitly Fig. [9] for se- 
lected optical indices. Both Fel 5270A and Mg2 present a mild cor- 
relation with [Fe/H] that, to a lesser degree, also applies to the G- 
band. The M05 models are sensitive to metallicity, and they do re- 



produce the observations. It is worth mentioning that there are two 
outliers in Fig. [8] NGC 6528 and NGC 6864. The former is located 
close to the Galactic centre and, thus, Wh/3 may be attenuated by 
reddening. Regarding NGC 6864, there is a s ignificant contributio n 
of blue horizontal branch stars (see Fig. 5 of lKravtsov et alfeOQTl) , 
which may be enhancing the H/3 absorption. 

The perception gained when considering absorption features is 
that M05 models, which up to date are the most suitable to describe 
almost all absorption lines/bands observed in the NIR, do properly 
reproduce the optical absorption line strengths. Regarding the NIR 
absorption lines, the models do underestimate the strengths of Mg I 
1.49/im, but they can properly reproduce the observed W\ of Fe I 
1.58^im, Si I 1.59^im, and CO 2.29/im, in about half of our sample. 

In the case of C O 2.29/im W\, similar results were obtained 
bv lLvub enova (2010, see their Fig. 9) for old and metal-poor Large 
Magellanic Cloud (LMC) globular clusters. These authors also 
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show that a large fraction of carbon stars would mimic the spec- 
trum of a younger cluster. To check this hypothesis, we overplot 
in Fig. [7] 1 Gyr models (asterisks joined by dot-dashed line). In- 
terestingly, the intermediate age models can describe the observed 
Wco 2.29pm in about half of our sample. Thus, we suggest that 
the NIR light in the core of these clusters (NGC 104, NGC362, 
NGC 2808, NGC 6388, NGC 6528, NGC 6553) may, perhaps, be 
dominated by C-rich stars. Note that, the metallicity distribution 
of the clusters is heterogeneous (see Tab. O, thus these large val- 
ues of Wx of the CO 2.29/im are probably not related to metal- 
licity. This suggests a possible carbon-fraction/age deg eneracy that 
affects NIR stellar populations. Indeed, according to Marig o et al.l 
(2009) (see their Fig. 5), the K luminosity fraction of TP-AGB 
stars in 10 Gyr old SSPs of different metallicities has been shown 
to grow from about 15% (at [Z/Z Q ]=0.0) to 50% (at [Z/Z ]=-2.3). 
We remind that the 1 Gyr models in Fig.|7]are not related to cluster 
age. Instead, thei simply reflect the presence of C-rich stars in some 
globular clusters. Also, it is interesting to remark that C-rich stars 
are only important in the NIR, beeing essentially unnoticed in the 
optical (see Fig.(8](. 

Interestingly, the Fel 1.58 fim atomic lines are well re- 
produced by the 1 Gyr population in about half of the sample 
(NGC 362, NGC 6388, NGC 6528, NGC 6553 and NGC 6864). In 
the case of Mg I 1.49/im, it occurs for almost all objects (NGC 104, 
NGC 2808, NGC 6388, NGC 6528 and NGC 6864). This suggests 
that the inclusion of empirical spectra of C- and O-rich stars in the 
models may also have an important effect in the strength of metal- 
lic lines. Another possibility lies in the fact that we are studying 
Q-enhanced Galactic globular clusters and, as the models do not 
include a-enhancement, this might explain why the models cannot 
reproduce the Wx of absorptions involving a-elements. 

Another powerful way to test SSPs models is by comparing 
observed and predicted colours. For this purpose, we collected the 
optical to NIR colours of our star clusters from the literature. These 
values are listed in Tab. [2] together with some basic properties of 
the cluster sample. 

We compare the observed colours with the models prediction 
in Fig. [10] Note that the optical colou r s wer e reddening-corrected 
using the ICardelli. Cla yton & Mathis (1989) extinction law, and 
the E(B-V) values from lBica et alj J2006h . Clearly if standard er- 
rors, in colours, are taken into account, the models are able to de- 
scribe almost all the observations. However, there are some outliers 
(NGC 6440, NGC 65 17, NGC 6528 and NGC 6553) for which it is 
not possible to determine whether or not the models apply, or if 
the photometric uncertainties are large, since these objects are lo- 
cated very close to the Galactic centre, and thus are more likely 
contaminated by field stars. Note that these objects are located in 
the Galaxy bulge (see Tab.O, and therefore, precise measurements 
are difficult. In addition, uncertainties in SSPs due to stellar tracks 
and EPS codes in t he V — K colour are typically 0.25 mag (see 
ICharlotetal.il 19961 M05). In this way, we suggest that M05 mod- 
els, in general, are able to predict the observed color values for the 
optical. For the NIR, however, they do systematically predict bluer 
colours, and a 2 Gyr SSP is required to properly reproduce the ob- 
served NIR colours. Thus, they may underestimate the age for the 
old population. 

The above results lead us to conclude that M05 models can 
provide reliable information on the NIR stellar population of galax- 
ies only when Wx and colours are taken together, in other words, 
W\ and continuu m fluxes should be simultaneously fitted (e.g 
iRiffel et all 120081) . alternatively the whole underlying spectrum 
should be used (e.g. IRiffel et~al]|2009l [201ol : iMartins et alj|2010h . 



However, the results should be taken with caution, since the mod- 
els tend to predict results biased towards young ages. A more robust 
test is fit the whole globular clusters spectra with models. However, 
this is out of the scope of this paper and we leave it for a forthcom- 
ing publication (Ruschel-Dutra et al., 2010 in preparation). 

4 CONCLUSIONS 

We use NIR SOAR/OSIRIS integrated spectra of 12 Galactic glob- 
ular clusters to test M05 NIR EPS models, and to provide spectral 
observational constraints to calibrate future NIR EPS models. M05 
models are used for being, up to date, the most suitable to describe 
absorption lines/bands in the NIR. Our main conclusions are: 

• The spectra of the cluster sample appear qualitatively similar 
in most of the NIR absorption features. 

• Many atomic absorption features like: A 1.49/im Mgl, A 
1.58/^m Fe i/Mgi, A 1.59/xm Si I, A l.llpm Mg I, A 2.21/xm Na 
I and A 2.26/im Cal as well as the A 1.62/xm, A 2.29/im CO and 
A 2.05/im CN molecular bands are clearly detected and identified 
in the spectra. The Wx of these features were measured, as well 
as the optical Wx of of G band (4300A), MgH (5102A), and Fel 
(4531 A). The globular clusters observations (W\ and colours) were 
compared with models predictions with ages from 4 up to 15 Gyr, 
with metallicities between ^jj ZQ and 2 ZQ. 

• M05 models are able to properly reproduce the optical Wx 
observed in globular clusters. 

• The models do underestimate the strength of Mg I 1 .49 pm, but 
they can properly reproduce the observed Wx of Fel 1.58/im, Si I 
1.59/im and CO 2.9/im in about half of our sample. For the remain- 
ing objects, we needed to consider intermediate-age populations. 
Thus, we suggest that the presence of C- and O-rich stars in the 
models is important to reproduce the observed strengths of metal- 
lic lines. Another possibility is the lack of a-elements enhancement 
in the models. 

• No significant differences are observed in the prediction of 
blue horizontal branch and standard models in the case of the Fe I 
lines (optical and NIR). A similar trend is observed for the Nal 

5 895 A. While in the case of the G-band the models which include 
blue horizontal branch do describe better the observations. No dif- 
ferences between blue horizontal branch and standard models are 
found for the NIR observations. 

• In general, M05 models can reproduce the observed colours in 
the optical, while in the NIR they tend to underestimate the age of 
the old population. 

The NIR spectral region is the most convenient one to study 
the stellar content of highly obscured sources. Thus, well calibrated 
and accurate NIR models is fundamental. Besides, observations us- 
ing adaptive optics, to properly correct for atmosphere effects, have 
become commonly used in th e NIR region, the reby allowing for 
high quality spectroscopy (e.g. IRiffel e t al. 2010). The availability 
of detailed, high spectral resolution, and well calibrated stellar pop- 
ulation models in the NIR will open a new window to disentangle 
the NIR stellar content in the inner few pc of nearby galaxies. Thus, 
the data presented in this paper will be important for the new gen- 
eration of NIR EPS models. 
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